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Abstract 

It has been known that Little Higgs models with T-parity, which can give a dark matter 
candidate, suffer from the anomalies of given models through the Wess-Zumino-Witten 
term, which in turn can violate the T-parity. Here we thus introduce a new kind of 
discrete symmetry called Z-parity, which is a residual gauge symmetry of an additional 
U(l) symmetry of the extended gauge group of the Little Higgs models. Because Z-parity 
comes from the gauge symmetry, this discrete symmetry not only remains unbroken under 
the anomalies but also can give the lightest Z-parity particle, which is stable and becomes 
a viable dark matter candidate. We also show that there is an allowed parameter space in 
the global symmetry breaking scale / versus the mass of the dark matter, which satisfies 
the current relic density of the cold dark matter. 
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1. Introduction 

At one time there had been a wisdom claiming that the supersymmetry is the only 
candidate to solve the hierarchy problem without big tuning since the supersymmetry 
itself gives the beautiful cancelations between boson and fermion loops and makes these 
nervous quantum corrections mild. According to this wisdom, the cancelation is based on 
the opposite sign between boson and fermion loops in Feynman rules. However, this kind 
of cancelation is not unique, that is to say, the possibility that particles with the same 
statistics can cancel each other's contribution under quantum corrections may exist. 
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Actually in the Higgs loops of the minimal supersymmetric standard model (MSSM) 
some diagrams among boson loops cancel each other's loop correction. 

Recently one novel and different approach for solving the hierarchy problem was sug- 
gested by Arkani-Hamed, Cohen and Georgi They revived the old idea that the Higgs 
boson is a pseudo-Goldstone boson from a spontaneous symmetry breaking of approxi- 
mate global symmetry Q], and made one successful model called the Littlest Higgs model, 
which can protect the electroweak scale under quantum loop corrections. In their first 
successful model, unlike supersymmetric cases, there are only a few additional new gauge 
bosons and a vector-like top quark to cancel their corresponding one-loop corrections of 
the partners with the same statistics. This way to cure the harmful quadratic divergences 
of the squared Higgs mass is surely more economical and simpler than that of the MSSM 
which has many more free parameters. Subsequently, many authors have been inspired 
by these attractive ideas, and have built various models [3J, which are largely classified 
by the breaking patterns of global symmetry. 

As soon as Little Higgs models had been introduced, these models also received a 
great attention from a different point of view besides the stability of the electroweak 
scale: Cheng and Low showed that these models can have a discrete symmetry called 
T-parity like the R-parity in the MSSM Q , and this T-parity can also give a stable dark 
matter candidate by assigning the Standard Model (SM) fields to even and heavy new 
particles to odd parity, respectively. Moreover, because the T-parity makes dangerous 
quantum corrections weak by prohibiting some vertices in Feynman rules, it can make 
the Little Higgs models more consistent with severe electroweak precision constraints 
than before. However, although the T-parity has good advantages, lately Hill and Hill 
5] showed that the T-parity can suffer from anomalies and be broken because of oddness 
of Wess-Zumino-Witten (WZW) term under T-parity. This term contains information 
of anomaly physics of a given theory in the Lagrangian. 

In this letter we thus introduce a new discrete Z2 symmetry called Z-paritjj^. This 
new parity generally comes from a residual discrete gauge symmetry of an additional 



U(l) gauge group 



ft 



In order to realize Z-parity we must simultaneously consider the 



In our first version of manuscript, we called it the X-parity. However, later we found that the term 
"X-parity" has already been used in Ref. [(J. We thanks the authors of [(J for letting us know. 



anomaly-free conditions for the extended gauge groups as well as for the relevant Yukawa 
terms. As a specific example, we choose the minimal Little Higgs model called Mini Little 
Higgs model (MLH) [8], only because of the simplicity of the gauge structure SU(2) L x 
U(l)j x U(l) 2 . We extend MLH and show how to realize Z-parity. We also show that 
Z-parity odd particles are right-handed Majorana particles, and there exist at least two 
Z-parity odd particles in the minimally extended model, therefore, the lightest Z-parity 
odd particle can be a viable dark matter candidate. It is also interesting that, as shown 
in Table 3, the charge assignments of Z-parity odd particles show Z2 symmetries between 
U(l) 1 and U(l) 2 , as well as among particles in each U(l) charge. Especially, these Z- 
parity odd particles can interact only with the Z' gauge boson and scalar particles. 

In addition, we investigate the allowed regions of the parameters which satisfy the 
current relic density of the cold dark matter, and find very narrow allowed regions. It 
is important that there exists an upper bound of the global symmetry breaking scale, 
/ < 1.1 TeV, depending on the dark matter mass (see Fig. 1). We also show that the 
mass of top partner if cannot be larger than around 3 TeV by using the constraint of / 
value and, therefore, the range of the mass of t 1 is about 1.4 — 3 TeV (see Fig. 2). 

2. Realization of Z-parity in Little Higgs Model 
A. Aspects of the Mini Little Higgs model 

MLH is based on the SU(2) L x U(l) 1 x U(l) 2 gauge symmetry, and an approximate 
U(3) global symmetry which is broken to U(2). It is important to notice that these U(l) 
charge assignments, which protect the squared mass of the Higgs particle from radiative 
quantum corrections of charged weak gauge bosons W ± and Z° in the SM, are necessary 
ingredients in MLH because of the absence of new charged particles. These heavy charged 
partners generally cancel the quantum corrections of the W ± , Z° in usual Little Higgs 
spirit. (Note that there is no additional SU(2) gauge group in MLH.) However, because of 
the simplicity of the gauge structure we can easily obtain simple anomaly-free constraint 
equations for each U(l) charge assignment and extend the model to have the new residual 
discrete gauge symmetry. We call this parity Z-parity from now on. 

Now we briefly present the lagrangian for the original MLH, 

4£5£ - \d^\ 2 , (i) 
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^Yukawa = Vi QlHIr + X b Q L Hb R + X e E L He R + \ v E L Hv R + h.c (2) 
+ yi$LSt R + y 2 fipLipR + h.c , 

where is a triplet of the global U(3) and decomposed into an SU(2)l doublet H and 
an SU(2)l singlet S as <j) T = (H, S) T . The second line describes the usual Yukawa terms 
plus the coupling term of the right-handed neutrino v R for the see-saw mechanism, and 
the third line presents the mixing between the top quark and heavy top partner ip. Note 
that the existence of the same coupling y\ for the first terms of the second and third 
lines, and the coupling 1/2 for the additional colored vector-like quark ipL.R, are necessary 
to cancel the dangerous top loop corrections to the Higgs h field. (See Ref. [8j for details 
of MLH.) 

Table 1: Examples of anomaly-free U(l) 1 and U(l) 2 charge assignments for li = —hi = — ^ (upper 
two rows) and li = hi = — ^ cases (bottom two rows) where I2 is a zi charge for El, and hi is a zi 
charge for H and we assume that si, which is a zi charge for singlet S, is — ^. In order to get 21 charge 
assignments we use hypercharge relation z\ + zi = Y. 
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With the help of those Yukawa terms, each of which give own constraint equation for 
each U(l) charge, and anomaly-free constraint equations from the triangular anomalies 
under SU(2) L x 11(1^ x U(l) 2 gauge group, we can obtain anomaly-free charge assign- 
ments. We present two examples, l 2 = — /12 = — \ case and I2 — h,2 — — \ case, in Table 
1. Here I2 is a Z2 charge for E L and /12 is a Z2 charge for H . The first case shows an 
anomaly- free set shown in Ref. [8[ . For the details of all Yukawa constraints and anomaly 
cancelation conditions, please see the Appendix of this article. 

B. Realization of Z-parity in the modified Mini Little Higgs model 

Let us start by introducing an additional singlet scalar S", which has the mass around 

the cutoff energy scale of the model, and two right-handed Majorana fermions {X\, X 2 ) 
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to achieve the new discrete symmetry as the residual gauge symmetry. These additional 
particles and charge assignments under gauge group SU(2) L x U(l) 1 x U(l) 2 are presented 
in Table 2. Note that we cannot yet fix each U(l) charge in Table 2 since there are 
many possible solution sets, which satisfy necessary constraint equations from Yukawa 
terms and anomaly-free conditions. We will show why we need at least two Z-parity 
odd particles X\, X 2 and two singlet particles S, 5", through details of these conditions 
below. 

Table 2: Additional fields and quantum numbers in the modified Mini Little Higgs model (mMLH) with 
Z-parity. Each zi[F] and Z2[F] denotes each U(l) charge for a field F. S' is a heavy scalar particle 
around the cutoff energy scale. X\ and Xi are right-handed Majorana fcrmions. Y is the Standard 
Model hypcrcharge and Q is the electromagnetic charge. 
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In the anomaly structure, X particles can contribute only to B 3 and B4 conditions 
(see below) because of darkness (neutral and not weak-charged) of X particles. Hence, 
we can easily show that with a single X particle one cannot satisfy these full anomaly-free 
conditions due to the complete separation between the anomaly-free MLH sector and the 
X sector. Therefore, at least two X particles are needed to cancel each other's anomaly. 
Therefore, we assume here that the X sector consists of the two right-handed Majorana 
particles, X\ and X 2 . Note that introducing of two extra right-handed Majorana fermions 
is the minimal way to realize Z-parity in this modified model. 

From the linear B 3 condition (of Table 6 in Appendix) each Xi must satisfy 

z 2 \X x \+z 2 \X 2 \ =0 . (3) 

Therefore, z 2 [JTi] = —z 2 [X 2 \. It is interesting that this charge assignment automatically 
satisfies other cubic anomaly condition B4, 

z^] 3 + z 2 [X 2 } 3 = . (4) 
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However, this charge assignment has inconsistency if, in addition, we consider mass terms 
of X\ and X 2 through the Yukawa couplings with only one singlet S, 

£ = 4£S£ + ^Yukawa + A*,iS (X^C" 1 ^) + \ X , 2 S {X? X 2 ) + h.c , (5) 

because each X has the same Yukawa constraint equation 

Y x : z 2 [S]+2z 2 [X h2 ] =0 . (6) 

Here each z 2 [Yi i2 ] must have the same charge, i.e. — z 2 [S]/2 — 5/6, which contradicts 
Eqs. (3,4). We can thus conclude that the modified Mini Little Higgs model (mMLH) 
with a single scalar S and Yukawa couplings cannot satisfy these anomaly free conditions 
simultaneously. Therefore, we have to introduce one more singlet particle S', so that each 
mass term for Xi and X 2 particles can be generated through separate Yukawa interactions 
of the two singlets S and S'. It is important to understand that these anomaly cancelation 
conditions and additional Yukawa terms are keys to assure a new discrete Z 2 symmetry, 
called Z-parity, as a residual gauge symmetry. With this simple prescription the Yukawa 
couplings are rewritten by 

£ = + ^Yukawa + *X,lS (X^C^Xj) + \ X , 2 S> (XfC^Xt) + h.C , (7) 

and Yukawa constraints are, 

Y x : z 2 [S]+2z 2 [X 1 ] =0 , z 2 [S'] + 2z 2 [X 2 ] = 0. (8) 

These equations show that if z 2 [S'\ = — z 2 [S] condition is satisfied, then all anomaly 
conditions are automatically fulfilled. Please note that these Yukawa terms can cause 
Higgs particle to have a large loop contribution from X\ particle via the coupling between 
singlet S and X\: The loop effect can revive the previous hierarchy problem. However, 
due to the j suppression factor from the nonlinear parametrization for singlet S, we 
can overcome this quadratic divergence for sufficiently large / values. We will discuss 
possible parameter regions for /, which satisfy the dark matter relic density, in Sec. 3. 

Our final results for mMLH are presented in Table 3. Note that in this table we 
re-scale the U(l) charges for particles by multiplying by factor |. It is now worthwhile 
to notice that these charge assignments show interesting Z 2 symmetry relations between 
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Table 3: Quantum numbers for the singlet scalar particles S and S' , and the right-handed Majorana 
fcrmions Xi and X2 in order to cancel the U(l) 3 and U(l) — [gravity] 2 anomaly conditions. Note that 
all U(l) charges are re-scaled by 4 to show simple charge relations between Xi, X2, S and S' . 
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U(l) 1 and U(l) 2 , as well as between X particles and between scalars S and S' for each 
U(l) charge. It is also remarkable that these particles can interact with Z' gauge bosons 
under the assignments of U(l) charges due to different signs of the charge values. 

Before we define Z-parity, we digress to consider the source of the extra singlet S'. 
This singlet S' may come from various origins. For example, this S' may be a complete 
singlet in low energy effective theory involving new physics at TeV order, or a singlet 
which breaks an extended gauge group at high energy beyond the cutoff scale of the 
model. We briefly assume that the scalar 5" particle comes from the extended gauge 
group at high energies. We would like to comment on two possible cases for the extended 
gauge group: (I) additional U(l) from the symmetry breaking of the grand unified theory 
or more enlarged gauge group, or even possibly (II) the gauge group including left-right 
symmetry gauge group SU(2) L x SU(2) R . (Note that in previous subsection we showed 
that there exists one interesting charge assignment corresponding to the I2 = hi = — j 
case in Table 1, where only right-handed particles couple to Z' gauge boson except for 
the scalar particles and top partner t' .) 

Finally, we present the long-awaited Z-parity from our anomaly-free charge assign- 
ments. In mMLH the singlet S breaks the original global U(3) into U(2). This singlet 
S also mixes U(l)i with U(l)2, so that we have the hypercharge gauge symmetry U(l)y 
and broken U(l)', which has the heavy Z' gauge boson. After U(l)' charges are nor- 
malized to integers by multiplying by a suitable constant number, we get the following 
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Tabic 4: Examples of U(l)' charge assignments and Z-parities (Z) for I2 = — h<z = — ^ (upper two rows) 
and h = h2 = —\ (lower two rows) cases. Note that all U(l)' charges are normalized to integers. 
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discrete Z 2 parity, which is denoted by Z, 

Z [all Fields except X] = even, Z [X] = odd. (9) 

Two examples are presented in Table 4. As can be seen easily, X particles are the only Z- 
parity odd particles under the discrete Z 2 symmetry. Therefore, one of these X particles 
can become a viable dark matter candidate. 

3. Dark matter candidate and the relic density 

Let us now consider the dark matter relic density for the lightest X particle. We have 
basically four free parameters, i.e. the masses of X particles (Mx x and Mx 2 ), the global 
symmetry breaking scale (/) and the Yukawa couplings (yi or y 2 ). Since X 2 particle 
is much heavier than X\ (note that X 2 particle has the mass via non-zero VEV of the 
5', which breaks the high energy gauge group), X\ can become the cold dark matter. 
We omit the mass parameter Mx 2 from the free parameter space because this heavy X 2 
particle does not contribute to the annihilation channel of relic density at the energy 
scale below a few TeV. We choose Mx i: f and y\ as free parameters and investigate 
allowed regions which can satisfy present relic density, in the interval 

550 GeV < M Xl < 2000 GcV, 500 GcV < f < 2000 GeV, 0.8 < yi < 4.0 

for fixed M X2 = 10 TeV. 
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The current dark matter abundance from WMAP collaboration Q is 
0.1037 < Q Xl h 2 {la) < 0.1161 (0.1099 ± 0.0062(lcr)) . 



To probe the allowed regions, we use the micrOMEGAs_2.2 [10J, which calculates the relic 
density as well as direct and indirect detection rates of the dark matter. First we show the 
allowed regions of the parameter space for / and Mx x in Fig. 1, which predicts / < 1.2 
TeV for the allowed range of Mx x < a few TeV. Here the red (light) and the blue (dark) 
regions represent allowed regions which satisfy the present relic density at 95% (2a) and 
99% (3a) confidence levels, respectively. It is interesting that these allowed regions look 
like step functions with step-like increase around 500 — 550 GeV and 900 — 1000 GeV. In 
fact, this behavior is closely related to the opening of new annihilation decay channels, 
X\Xi —¥ hS and X±Xi — > SS. In the first step, the opening of the new annihilation 
channel, hS, entirely dominates annihilation of the dark matter. Namely, as the relative 
contribution of the new annihilation channel increases, the relative contribution of the 
previously dominant W + W~ and ZZ decay channels decrease. Similarly, in the second 
step new annihilation decay channel SS begins to open when the mass of dark matter 
particle Xi approaches about 880 GeV. As the contribution of the new SS decay channel 
gets larger, the relative contribution of the previous dominant decay channel hS rapidly 
decreases. Eventually the contribution of the SS decay channel totally governs decay 
modes of the dark matter. 

Because the top quark contribution in the relic density is so negligible, we may expect 
that ?/i parameter can also be omitted from the four free parameters. However, the value 
of j/i is important to relate the mass of top with top partner particles. From the Yukawa 
terms (See Eq. (2)) we can calculate the relevant masses by 



where mt is mass of the top quark and m t < is mass of the top partner. In Fig. 2 we show 

the lower bound of the mass of top partner t' as a function of symmetry breaking scale 

/. It is interesting that the mass range of the top partner t' is about 1.4 — 3 TeV for the 

allowed values of / < 1.2 TeV, which is estimated from Fig. 1. 

We note that recently there have been interesting observations from PAMELA [ll| 

and ATIC [12j experiments, which support a very striking energy source of energetic 
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Figure 1: The allowed regions which satisfy the current relic density of the dark matter X\. The blue 
(dark) and red (light) parts correspond to 99% and 95% confidence levels, respectively. The / is a global 
symmetry breaking scale and Mx j is mass of dark matter particle X\ . Note that the allowed region is 
very narrow and there are also upper and lower limits of / along Mx-i ■ 
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electrons and positrons within ~1 kpc away from the Solar System. At present this 
source has not yet been fully identified. The source may be the annihilation of dark 
matter particles or an astrophysical object. Some papers show that these objects may be 
pulsars which produce a bump in the spectrum [13j, and recent results from the Fermi 
Telescope [l4| suggest that the ATIC result may not be a signal of the annihilation 
of the dark matter particles. In mMLH with Z-parity, these X particles seem to be 
difficult to have selective annihilations to charged leptons in order to produce the excess 
of electrons and positrons. Consequently, the X particles are not the sources of the ATIC 
and PAMELA observations. 

Finally, we remark that the modified Mini Little Higgs Model has two interesting 
implications: 

(I) Various Little Higgs models have extended gauge symmetries, e.g., some models 
have the extended hypercharge gauge symmetry like U(l) 1 x U(l) 2 before the global 
symmetry breaking. Therefore, we can apply our Z-parity to these kinds of general 
Little Higgs models, and we can infer that there can exist stable particles other than 
usual T-parity stable particles. 

(II) Another interesting aspect is a possibility to explain the neutrino masses, cold 
dark matter (CDM) and hierarchy problems simultaneously. There are papers which 
explain the CDM and neutrino mass simultaneously but not the hierarchy problem [15I ] . 
However, mMLH can not achieve above aspect because of the minimal structure of the 
scalar sector of MLH. 

In conclusion, we showed how to realize Z-parity in the modified Mini Little Higgs 
model, and found that Z-parity can give a new kind of stable and neutral dark mat- 
ter candidate, the right-handed Majorana particle. Some additional fields S' and X2 
particles were also introduced to accomplish the anomaly-free conditions and to obtain 
simultaneously the Yukawa terms leading to a residual gauge symmetry in the model. 
Moreover, we also found that the dark matter X\ can satisfy the present dark matter relic 
density, and there exist allowed regions in the free parameter set of the model. Finally, 
the mass range of the top partner t' was also deduced for the allowed values of the global 
symmetry breaking scale /. 



11 



ACKNOWLEDGEMENTS 

C.S.K. was supported in part by Basic Science Research Program through the NRF 
of Korea funded by MOEST (2009-0088395) and in part by KOSEF through the Joint 
Research Program (F01-2009-000-10031-0). JBP thanks S. K. Kang for valuable discus- 
sions and encouragement about this paper and H. S. Lee for his kind and critical opinions. 

APPENDIX. 

In this appendix, we briefly introduce the particles and corresponding quantum num- 
bers of MLH, and then discuss the Yukawa terms and anomalies in order to help under- 
stand of the main body. Subsequently, we show how to obtain the anomaly-free charge 
assignments in MLH, and present anomaly-free solution sets. 



Table .5: Fields and quantum numbers in the Mini Little Higgs model. Each Zi[F] denotes U(l) charge 
for a field F. Here tpL.R is a vector-like colored top partner, and 5 is a scalar particle that breaks the 
global U(3) symmetry. Y is the Standard Model hypcrcharge and Q is the electromagnetic charge. 
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First, the particles and corresponding quantum numbers under SU(3) C x SU(2) L x 
U(l) 1 x U(l) 2 are presented in Table 5. The first column is the symbol of particle. The 
second and third columns are the SU(3) C and SU(2) L charges for color and chiral weak 
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interaction, respectively. Next two columns are the U(l) charges of the corresponding 
field denoted by Zi[F], where i is the index of each U(l) gauge group. Note that we 
have not yet assigned the U(l)i and U(l)2 charges because we have to discuss various 
constraints from Yukawa terms and anomalies of the given model. These U(l) charge 
assignments protect MLH from the radiative quantum corrections up to 5 — 10 TeV scale. 



Table .6: Anomalies and anomaly cancelation conditions. Here Aj's denote the anomalies involving 
U(l)j only. Similarly, Bj's denote the anomalies involving U(l) 2 only. Cj's are used as the mixed 
anomalies involving U(l)j and U(l) 2 . 
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We summarize all important anomalies in Table 6, and divide all these anomalies into 
three classes (A, B and C). Here Aj's are the anomalies involving U(l) 1 only. Similarly 
Bj's denote the anomalies involving U(l) 2 only. The mixed anomalies between these 
U(l) 1 and U(l) 2 are denoted by Cj's. 

Now we consider the solution sets of the U(l) charges. In MLH we have 18 unknown 

parameters such as z i= i^ 2 [F], where F — {Ql, ta, &r, El, vr, en, tpL,R, H, S}. There are 

also 10 = 2 x 5 constraints : 5 Yukawa constraints (Yu,Yd,Ye,Yn,Ys) per each U(l), 

where subscripts denote up, down, electron, right-handed neutrino and the singlet S, 

respectively. Moreover, there are 9 hypercharge relations and 10 anomaly conditions 

(Ai_4, £?i_4, C1.2). Therefore, the number of constraints is larger than that of unknown 
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parameters. However, all these constraints are not completely independent. For instance, 
anomaly conditions ^2,3 and i?2,3 give the constraint equations which are the same as 
Yukawa constraints given above. In order to find the solutions of the valid charges 
for U(l) 1 and U(l) 2 efficiently, it would be better to use hypercharge relations first. 
These 9 hypercharge relations reduce our original 18 unknown parameters to 9 unknown 
parameters, and simplify two (Ci and C2) anomaly cancelation conditions by adding and 
subtracting these two conditions. Then, we can show that these two nonlinear equations 
become simpler, 

/ 

'2. _„ ( 2 ) 



Ci 
C2 



/ 

Since 5 Yukawa constraints and B\ give a total of 6 linear constraints for 9 unknown 
parameters for the U(l) 2 sector, three parameters still remain to be unknown. Therefore, 
general solution sets of U(l) 2 are given by 



Z2 [Q L ] 
Z2 [*«] 
Z2[bR] 
Z2[E L ] 
Z2[vr] 

Z2 [en] 

z 2[lpL,R] 
Z2[H] 

\ z 2 [S] ) 



h 
3 



f-l\ 
-1 
-1 
3 
3 
3 

-1 


V / 



/ 



\ 
1 

-1 

1 

-1 
1 
1 

V / 



*2 








1 




(.3) 



where I2 = 22 [El], ^2 = Z 2[H] and S2 = 2^2 [S]- It is interesting that these general solu- 
tion sets automatically satisfy remaining constraints, quadratic Ci, C2 and cubic B4. We 
present here two example assignments for I2 = h% = \ and I2 = — 112 = 9 cases in Table 1. 



[1] N. Arkani-Hamed, A. G. Cohen and H. Georgi, Phys. Lett. B 513, 232 (2001) 
[arXiv:hep-p h/0105239 ; N. Arkani-Hamed, A. G. Cohen, T. Gregoire and J. G. Wacker, JHEP 
0208, 020 (2002) [arXiv:hep- ph/0202089 ; For a review and references: T. Han, H. E. Logan, 
B. McElrath and L. T. Wang, Phys. Rev. D 67, 095004 (2003) |arXiv:hep-ph/0301040] ; M. Schmaltz 

14 



and D. Tucker-Smith, Ann. Rev. Nucl. Part. Sci. 55, 229 (2005) [arXiv:hep- ph/0502182 ; 
M. C. Chen, Mod. Phys. Lett. A 21, 621 (2006) |arXiv:hep-ph/0601126] . 
[2] H. Georgi and A. Pais, Phys. Rev. D 10, 539 (1974); H. Georgi and A. Pais, Phys. Rev. D 12, 508 
(1975). 

[3] N. Arkani-Hamed, A. G. Cohen, E. Katz and A. E. Nelson, JHEP 0207, 034 (2002) 
arXiv:hep-ph/0206021 ; N. Arkani-Hamed, A. G. Cohen, E. Katz, A. E. Nelson, T. Gregoire 
and J. G. Wacker, JHEP 0208, 021 (2002) [arXiv:hep-ph/0206020j ; W. Skiba and J. Tern- 
ing, Phys. Rev. D 68, 075001 (2003) [arXiv:hep-ph/0305302| ; D. E. Kaplan and M. Schmaltz, 
JHEP 0310, 039 (2003) [arXiv:hep-ph/0302049| ; M. Schmaltz, JHEP 0408, 056 (2004) 
[arXiv:hep-ph/0407143' ; E. Katz, J. y. Lee, A. E. Nelson and D. G. E. Walker, JHEP 0510, 088 
(2005) [arXiv:hep-ph/0312287| ; K. Agashe, R. Contino and A. Pomarol, Nucl. Phys. B 719, 165 
(2005) [arXiv:hep-ph/0412089] . 

[4] H. C. Cheng and I. Low, JHEP 0309, 051 (2003) [arXiv:hep-ph/0308199| ; H. C. Cheng and 
I. Low, JHEP 0408, 061 (2004) |arXiv:hep-ph/0405243| ; I. Low, JHEP 0410, 067 (2004) 
[arXiv: hep-ph/0409025] . 

[5] C. T. Hill and R. J. Hill, Phys. Rev. D 76, 115014 (2007) |arXiv:0705!0697l [hep-phi 1 . 

[6] A. Freitas, P. Schwaller and D. Wyler, JHEP 0912, 027 (2009) |arXiv:0906.T8T6l [hep-ph]] . 

[7] L. E. Ibanez and G. G. Ross, Nucl. Phys. B 368, 3 (1992); R. N. Mohapatra and M. Ratz, Phys. Rev. 
D 76, 095003 (2007) |arXiv:0707!4070l [hep-ph11 ; H. S. Lee, C. Luhn and K. T. Matchev, JHEP 0807, 
065 (2008) |arXiv:0712."3505l [hep-phll; H. S. Lee, Phys. Lett. B 663, 255 (2008) |arXiv :0802.0506 
[hep-ph]]; T. Hur, H. S. Lee and C. Luhn, JHEP 0901, 081 (2009) | arXiv:0811."08T2l [hep-ph]] . 

[8] Y. Bai, Phys. Lett. B 666, 332 (2008) | arXiv:0801.T662l [hep-ph[[. 

[9] J. Dunkley et al. [WMAP Collaboration], Astrophys. J. Suppl. 180, 306 (2009) |arXiv:0803!0586l 
[astro-ph]] . 

[10] G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 180, 747 (2009) 
arXiv:0803.2360 [hep-ph]]; G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. 
Commun. 176, 367 (2007) [arXiv:hep-ph/0607059] . 
[11] O. Adriani et al. [PAMELA Collaboration], Nature 458, 607 (2009) | arXiv:0810.4995l [astro-phll . 
[12] J. Chang et al. [ATIC Collaboration], Nature 456, 362 (2008). 

[13] S. Profumo, larXiv:0812.4457l [astro-ph], D. Hooper, P. Blasi and P. D. Serpico, JCAP 0901, 025 
(2009) |arXiv:0810T527l [astro-ph]]. 

[14] A. A. Abdo et al. [The Fermi LAT Collaboration], Phys. Rev. Lett. 102, 181101 (2009) 
arXiv:0905.0025 [astro-ph.HE]] . 

[15] J. Kubo and D. Suematsu, Phys. Lett. B 643, 336 (2006) [arXiv:hep-ph/0610006] ; T. Kajino, 
M. Kusakabe, K. Kojima, T. Yoshida, D. G. Yamazaki, K. Ichiki and G. J. Mathews, Mod. Phys. 
Lett. A 23, 2427 (2008); A. Anisimov and P. Di Bari, larXiv:0812. 50851 [hep-ph], H. Sung Cheon, 
S. K. Kang and C. S. Kim, JCAP 0805, 004 (2008) |arXiv:0710!2416l [hep-phll. H. Sung Cheon, 
S. K. Kang and C. S. Kim, larXiv:0807.0981l [hep-phi. H. S. Cheon, S. K. Kang and C. S. Kim, 
larXiv:0909.2T30l [hep-ph] . 



15 



